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ABSTRACT 

This paper describes the development of a mathe- 
matical model for predicting the strength and micro- 
mechanical failure characteristics of continuously 
reinforced ceramic matrix composites. The local-global 
analysis models the vicinity of a propagating crack tip 
as a local heterogeneous region (LHR) consisting of 
spring-like representation of the matrix, fibers and 
interfaces. This region is embedded in an anisotropic 
continuum (representing the bulk composite) which is 
modeled by conventional finite elements. 

Parametric studies are conducted to investigate 
the effects of LHR size, component properties, inter- 
face conditions, etc. on the strength and sequence of 
the failure processes in the unidirectional composite 
system. The results are compared with those predicted 
by the models developed by Marshall et al . (1985) and 
by Budiansky et al . (1986). 

INTRODUCTION 

The failure characteristics of fiber reinforced 
composites are dictated by various micromechanical 
failure processes such as matrix microcracking, slip- 
ping between matrix and fibers, delamination and fiber 
breakage. This paper presents a local-global model (it 
combines micromechanical and macromechanical analyses) 
which considers the vicinity of a crack tip a "process 
zone" capable of modeling such phenomena. Of special 
interest to our study are quantities such as critical 
matrix cracking stresses, since such damage leads to 
oxidation and eventually to fiber degradation, and is 
therefore used as an important criteria in design. 

Also of interest is the global response of the system 
to external loads (axial P-6 curves) and the ultimate 
load that the composite can sustain. 

Before presenting a detailed description of the 
present model, a comparative review of existing models 
for fracture mechanics of brittle matrix composites (BMC) 
is warranted. The most quoted models are those developed 
by Aveston et al. (1971), Marshall et al . (1985) and 
Budiansky et al . (1986). The following are some of the 
key concepts underlined in these models which we shall 
seek to study or validate through our model. 


Marshall et al . (1985) have developed a model 
which can be used to predict the stress at which a 
matrix crack propagates across the specimen. This 
stress will, henceforth, be called the critical matrix 
cracking stress, a m at- In their analysis, a frictional 
bond between fibers and matrix was assumed whereby 
slipping takes place when the interface shear stress 
reaches a critical value. Using a stress intensity 
factor approach, they have shown that a distinction 
needs to be made between short and long cracks. Short 
cracks are those for which the entire crack length con- 
tributes to the stress intensity factor as a result of 
fiber bridging, and therefore propagate at a stress 
which depends on the crack length. Long cracks experi- 
ence a crack mouth displacement which asymptotically 
approaches a constant value u 0 . This limiting dis- 
placement is reached at a distance c 0 from the crack 
tip. For such cracks, c ma t is independent of the 
crack length, since the contribution to the stress 
intensity factor from the fibers is limited to the 
length c 0 behind the crack tip. It is important to 
note that this model implicitly assumes the stress- 
strain diagram shown in Fig. 1(a), since no nonlineari- 
ties are assumed prior to the matrix cracking stresses. 
It will be shown using the model proposed in this paper 
that this assumption leads to a good estimate of a ma t- 
However, for various constituent properties the present 
model shows that significant nonlinearities may occur 
prior to ° m t ( Fig. Kb)). These irreversible defor- 
mations, which are due to slipping between fiber and 
matrix and microcracking may prove to be significant 
for fatigue types of loading. 

Budiansky et al . (1986) have considered steady 
state matrix cracking stresses for two conditions: 

(1) unbonded, frictionally constrained fibers, where 
the frictional restraint is the same as in Marshall 
et al. (1985); and (2) initially bonded fibers which 
debond due to crack tip stresses. The analysis is 
based on the Griffith energy criterion which considers 
the change in potential energy with respect to crack 
growth. The critical cracking condition is associated 
with the upstream and downstream stress states, far 
ahead of and behind the crack front. 
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For case (1), the results generalize those of the 
ACK theory by considering matrix cracking stresses for 
conditions which lie between the no-slip and the large 
slip cases. They showed that the critical cracking 
stress, a ma t» can be obtained using the graph shown in 
Fig. 2 in conjunction with Eqs. <1) and (2). The pro- 
cedure includes first evaluating the two parameters a o 
and a ] given by 


!o 

E = B 


6 »f E f 


v~;E( 1 + v ) 
m m 



G 


m 


rE 


L mj 


1/2 


(1) 


" Ev f E f T 0 

i / j 

M 

v EE 


rE 

m m 


L mj 


1/3 


where B is a utility constant given by 

o -,1/4 


B * 


2 v: 


-6 In (v.) - 3v ( 3 - v f > 
r m t 


( 2 ) 


Ef and E m are the Young's moduli of the fiber 
and matrix, E is the composite modulus using the rule 
of mixtures, v f and v m are fiber and matrix volume 
fractions respectively, r is the fiber radius, G m is 
the critical strain energy release rate of the matrix, 
and x o is the interface shear strength. 

It is observed that for a ]/ a o < 1.0, the ACK 
expression for large slip cracking stress is recovered, 
i .e. , 


THE LOCAL MODEL 

The local heterogeneous region (LHR) shown in 
Fig. 3 consists of three distinct components: fibers, 

matrix, and fiber-matrix interface. Following Kanninen 
et al . (1977), each component is assumed homogeneous 
and isotropic and is modeled by a succession of rectan- 
gular spring elements (Fig. 4). 

Each LHR spring element consists of four exten- 
sional springs and one rotary spring at each of the 
corner nodes. Each node has 2 degrees of freedom as 
shown in Fig. 5. An assembly of such elements behaves 
like a homogeneous anisotropic continuum in a state of 
plane deformation. 

The stiffness matrix of a LHR element is given by 

{F} = [ K ] { d } , (5) 

where {F} is the nodal force vector, [K] is the ele- 
ment stiffness matrix, and {d} is the nodal displace- 
ment vector. 

The LHR element spring stiffnesses are related to 
the elastic properties of the material. The deriva- 
tions for the spring stiffnesses and for the LHR stiff- 
ness matrix are given in Appendix A. 

Each LHR fiber and matrix element is capable of 
fracturing in either of the four possible modes shown 
in Fig. 6. Modes 1 and 2 correspond to crack growth 
in the x-direction by an amount a/2, while modes 3 
and 4 represent cracking in the y-direction within 
the element by a length of b/2, where a and b are 
the length and width, respectively, of a LHR element. 
These damage lengths have been used for the sake of 
simplicity in representing fracture events in the LHR 
elements. Each event contributes to a loss of stiff- 
ness of the element, and consequently leads to a change 
in the strain energy, AE S , given by: 
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where (^ iq ) is the Fracture toughness of the matrix. 

The ratio °] /°q is then calculated and Fig. 2 is used 
to obtain the corresponding value of a mat /a 0* 

For bonded-debonding fibers, the following expres- 
sion was derived: 
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where {u} is the displacement vector of the element, 
[SKI is the change in the element stiffness due to a 
fracture event, and Ac is the length of the crack in 
the element, i.e., a/2 or b/2. The above relation- 
ship is derived in Appendix B. 

It is assumed that for each material (fiber and 
matrix) a critical rupture energy is known. This pro- 
vides a decision rule for breakage in each separate 
element of the LHR. 

For purposes of comparison with the results of 
Marshall et al . ( 1985), preliminary simulations were 
conducted assuming that fiber slipping occurs at a 
critical shear stress In subsequent studies, more 

realistic modeling which includes delamination and an 
elastic-plastic friction model will be employed. 

THE LOCAL-GLOBAL MODEL 


where lq is the debond length and Gq is the criti- 
cal energy release rate of the interface. 

Expressions for the slip length l s (case <D) 
and debond length lq (case (2)) will be discussed 
subsequently. 

This paper presents the preliminary results 
obtained using a model which can be used to predict 
the fracture characteristics of reinforced ceramic 
matrix composites. This model has been used to simu- 
late an experiment in which a notched specimen is sub- 
jected to tensile stresses. The results are compared 
to those obtained using the aforementioned models. 


A schematic representation of a typical local- 
global model is shown in Fig. 7 (not to scale). The 
notched composite specimen modeled has a length 
1 = 20 cm and width w = 5 cm. The LHR surrounds the 
crack tip and is embedded in the bulk anisotropic com- 
posite, discretized by standard four-node constant 
strain isoparametric finite elements. The LHR and the 
outer zone are coupled through the enforcement of dis- 
placement compatibility at the nodes between the LHR 
and standard finite elements. 
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The specimen is loaded incrementally In tension 
perpendicular to the crack plane and parallel to the 
direction of the fibers. Displacements are evaluated 
at every nodal point in the LHR, and the strain energy 
associated with each possible rupture event in every 
LHR element is calculated. The critical rupture ener- 
gies of the fiber and matrix elements are obtained from 
their fracture toughnesses as follows 



where K 
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is the fracture toughness of the 


matrix (fiber) material, E 
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is the critical 


energy for a given length of damage growth within the 
matrix (fiber) LHR element; and v m (vf) is the Pois- 
son's ratio of the matrix (fiber). 

The critical regions, if any exist, are allowed to 
fracture in any one of the four modes described earlier 
and appropriate modifications are made in the LHR 
stiffness matrix. 

As stated earlier, failure of the interface ele- 
ments is based on the interfacial shear stress crite- 


ria. For each load level, the shear stress T 


xy 


is 


computed in every LHR interface element and this value 
is compared to the prescribed critical shear stress 


x q . If T X y < t q (no-slip condition) relative dis- 
placement between the fiber and matrix is constrained. 
For elements where T X y > x o> the fiber is allowed to 
slip through the matrix and appropriate changes are 
made in the shear stiffness of the interface element 
so that slip takes place at a constant shear stress 


T xy s T 0* 

The solution procedure is repeated for the next 
load increment to reflect additional local rupture 
events. The progression of the main crack through the 
matrix (and eventually through the fibers) is followed 
by expanding the LHR zone continuously. The load at 
which the LHR elements begin to fail without any fur- 
ther increment of load is recorded as the ultimate 
strength of the composite. This is generally found to 
occur at a stage when the main crack has broken through 
four (or sometimes less) successive fibers. 


RESULTS 


Physical Properties and Parameters 

Most of the results presented in this paper assume 
material properties (refer to Table I) obtained from 
the following sources: Weeton et al . (1986), DiCarlo 

(1984), and Bubsey et al . (1983). 

Two fiber volume fractions were considered, 0.2 
and 0.4. Consistent with these values were the proper- 
ties of the anisotropic bulk composite: 

v f * 0.2: Eyy = 242.8 GPa; E xx * 227.47 GPa; 

v X y = 0.3; G X y = 100.42 GPa 

v f = 0.4: Eyy * 279.6 GPa; E xx - 253.92 GPa; 

v X y = 0.3; G X y * 112.93 GPa 

The interface shear strength, t q, was varied in the 
practical range of 0 to 10 MPa. The fracture toughness 
of the matrix was also varied from 1.0 to 5.0 MPaV^* 


The diameter of the SiC fibers (2r) was taken as 
100 pm. Assumption of regular hexagonal packing arrange- 
ment gave an interfiber distance of about 213 and 150 pm 
(center to center) for volume fractions of 0.2 and 0.4, 
respectively. The interface was assumed to have a nom- 
inal thickness of 5 pm. 

To investigate the convergence of the local-global 
model, ultimate and matrix cracking stresses were stud- 
ied for various LHR sizes. Figure 8 shows the effect 
of starting LHR size on the ultimate stress a u it and 
the critical matrix cracking stress, a ma t* Convergence 
is observed for LHR lengths greater than about 1000 pm 
and hence LHR sizes of the order of 1000 to 1500 pm 
were used throughout the study. 

Comparisons With Results Obtained by Marshall et al . 
(1985) 


Effect of crack length . As mentioned previously, 
Marshall et al . (1985) have shown that when the crack 
length is greater than a certain characteristic size, 
matrix cracking stresses are independent of the crack 
length c. To check this result, critical matrix 
cracking stresses ( a ma t ) were calculated for various 

M 

crack sizes. For the case t q * 10 MPa, Kj C » 4.6 MPa^ 

and v f * 0.4, the results are presented in Fig. 9. 
Though” is high for small crack lengths, it is 

found to be independent of crack length for 
c > 0.046 cm. This compares favorably to the steady 
state crack length proposed by Marshall et al . (1985) 
obtained from the relation 
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M 

where n * VfEf/v m E m , is the fracture toughness of 

the matrix, and I is 1.2 for straight cracks. 

It should be noted that Eq. (8) was derived by M 
assuming that the matrix stress intensity factor, 1C , 

I i 

is related to the composite stress intensity, Kj , by 

the relation K* = K^E/E. This relation is not valid 
I I m 

unless the crack length is of the order of several 
fiber spacings. In a future communication, the range 
of validity of this equation will be reported. 

Effect of matrix fracture toughness on a mat and 

o ult . Having obtained the steady state crack length, 

the critical matrix cracking stress, 0 ma t> and ultimate 
strength of the composite, a u )t* were investigated next 
using an initial crack length such that steady state 
conditions are obtained. The ultimate strength and 
o ma t depend on several constituent material properties 
such as fracture toughness of the components, interface 
shear strength, fiber spacing and volume fraction of 
the fibers etc. Each of these parameters are consid- 

ered separately. The variation of a mat with K^ c 

over the range 1.0 to 5.0 MPa\/m is studied for a notch 
length of 1 cm and interfacial shear strength of 10 MPa. 
The critical matrix cracking stress is found to increase 
M 

with and the results compare well with those of 

Marshall et al . (1985) given by 
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The results from our model and those obtained from 
the above relationship are shown in Fig. 10. Since the 
results of Marshall et al. (1985) and Budiansky et al. 
(1986) are almost equivalent, they do not plot as sepa- 
rate curves. The ultimate strength is not affected 

significantly by variation in , since the ultimate 

strength of the composite is controlled primarily by 
the fiber-bundle strength. If ly is the length of 
the initially uncracked ligament in the notched speci- 
men of width w, the bundle strength is given by 
a f v f(l u /w), where °f is the tensile strength of the 
fibers. For the SiC fibers, a f = 1.83 GPa. Hence, 
for a notch length of 1 cm and w * 5 cm, the theoret- 
ical bundle strength is 586 MPa for v f = 0.4, and 
293 MPa for v f = 0.2. The ultimate composite strength, 
a ult’ obtained using our model is slightly less than 
these calculated values. This is due to the sequential 
breakage of the fibers leading to catastrophic failure. 

Effect of interfacial shear strength . Since 
interface shear strengths are difficult to estimate, a 
parametric study was conducted for t q ranging from 
5.0 to 10.0 MPa. The variation of a ma t with x q, 
shown in Fig. 11, shows good agreement with results 
obtained by Marshall et al . (1985) using Eq. (9). 


Effect of fiber volume fraction . Two fiber volume 
fractions vf were considered, 0.2 and 0.4. The fiber 
volume fraction controls the interfiber spacing and its 


M 


and 


effects on a mat and o ult with respect to 

T 0 are presented in Figs. 12(a) and (b). It is 
observed that a ma t and a u it values for v f = 0.2 
are approximately half of those for v f = 0.4. These 
results agree with those obtained using £q. (9). 


Variation of crack mouth opening displacement with 
crack length . Another key concept introduced in the 
model developed by Marshall et al. (1985), is that of 
the equilibrium crack opening u « u 0 for long cracks. 
To study the effect of crack length on u, the crack 
opening displacement at the mouth of the crack was 
recorded for various crack lengths (Fig. 13). Although 
u increases with increasing crack length for crack 
lengths less than 1 cm, it appears to approach a con- 
stant value of approximately 1 pm for crack lengths 
greater than that size. This value may be compared to 
the limiting displacement u 0 estimated by Marshall 
et al. (1985) by the relationship 

2 

u = 7 r (10) 

{ 4 Vf E f (1 ♦ n >} 

where a is the farfield applied load. For ° = a ma t 
* 300 MPa, x 0 = 10 MPa, and v f = 0.4, u 0 * 0.8 pm. 

The agreement is quite good. 

Comparison With Results Obtained by Budiansky et al . 

(1986) 


(Eqs. (1) to (3)). Because their results lead to pre- 
dictions equivalent to those of Marshall et al . (1985) 
they plot as the same curve in Figs. 10 and 11. 

Comparison of slip lengths . For unbonded compos- 
ites, the slip length has been derived by Budiansky 
et al . (1986) as 
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where l s is the slip length on either side of the 
crack face, a is the applied stress, and 
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For * 10 MPa, Kj C * 4.6 MPa^/m; and a = a ma ^ 

300 MPa, 1 * 806 pm. 

Results of our simulation are schematically pre- 
sented in Fig. 14. The slip length is approximately 
900 to 1000 pm on either side of the crack face, which 
agrees quite well with the result of Budiansky et al . 
(1986). 


Composite Failure Sequence 

While the existing models are concerned only with 
the steady state condition, the present model enables 
a detailed study of the failure sequence as the speci- 
men is loaded incrementally. A schematic representa- 
tion of the failure events in the LHR for a notch 

length of 1 cm, t q = 10 MPa, = 4.6 MPa \/m is shown 

in Fig. 14. Slipping first occurs at an applied stress 
of 140 MPa, and the interface continues to slip with 
increasing load. The first matrix crack is observed 
around the crack tip at a stress of 295 MPa. For 

matrices with lower K 1 ^, some microcracking is observed 

around the crack tip before the growth of the main 
crack. The matrix crack steadily progresses through 
the composite section and traverses it completely at a 
critical matrix cracking stress of about 305 MPa 
( a ma t). As loaci 1s increased, several secondary 
matrix cracks appear in the matrix. The first fiber 
failure occurs at a load of 505 MPa, and in general, 
this occurs at a load of about 90 percent of the ulti- 
mate stress. With increasing load, the crack breaks 
through three successive fibers before catastrophic 
failure occurs at an ultimate stress of 550 MPa. 


Load-Deflection Behavior 

The load-deflection curves as functions of the 
interfacial shear strength are shown in Fig. 15. 

Except for the case t q - 0, it is observed that non- 
linearities begin at point (a) before a ma ^ is 
reached. This decrease in stiffness is a result of 
fiber slipping and matrix microcracking. A further 
reduction in stiffness occurs at region (b). This 
region corresponds to extensive matrix cracking through 
the section of the composite, so that the load is now 
essentially carried by the fibers. The ultimate fail- 
ure load ( a u i t ) is reached at (c). 

CONCLUSION 


Comparison of critical matrix cracking stress 
a mat* The critical matrix cracking stresses obtained 
usTng our model are compared to those obtained by 
Budiansky et al . (1986) for unbonded composites 


A model is presented which can be used to predict 
the failure characteristics of fiber reinforced compos- 
ites. Preliminary results obtained using the model 
compare favorably with those predicted by existing 
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models. This model may be generalized to consider more 
complicated geometries and loading conditions as well 
as the behavior of composites containing a random dis- 
tribution of microflaws and strengths. It has been 
suggested (Dollar and Steif, 1988) that assumption of 
constant interface shear strength may not be accurate. 
The model is currently being improved through the 
implementation of more realistic interfacial constitu- 
tive modeling. In subsequent work, special friction 
interface elements developed by Plesha et al . (1987) 
will be used. 


APPENDIX A - DERIVATION OF LHR ELEMENT STIFFNESS 

For a homogeneous isotropic material, the stress 
strain relationship is given by 

R E R 

Vyy ) = (1 + v)(l - 2v) v 

y ^ ° 


v 0 

(1 - v) 0 

0 (1 - 2v) 
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The continuum is modeled by a set of spring ele- 
ments as described previously. The values of the 
spring constants are related to the material's elastic 
properties by the following relationships: 
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for all i , j = 


where a and b are the length (x-di rect ion) and 
width (y-di rection) of the LHR element respectively; 

K^ is the extensional stiffness in the x-direction 
between the i th and nodes; K^ is the exten- 

sional stiffness in the y-direction between the i th 

and j** 1 nodes; K^ is the cross extensional compo- 
xy 

nent (Poisson's contraction effect); and C. is the 
rotational stiffness at each node. 1 

The above relationships are derived as follows. 
Consider a finite continuum element of dimensions 
a x b, which is modeled by a LHR spring element of the 
same dimensions. By forcing the spring element to 
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behave like the continuum element, appropriate relation- 
ships may be obtained for the spring stiffnesses, K xx , 
K X y, Kyy, and Cj . Since the material is isotropic, 


.12 ^34 _ . . 

xx xx xx’ 


.14 .23 . 

K yy * K yy = K 


yy> 


element may now be formulated. The stiffness coeffici- 
ents can be evaluated through energy considerations as 
follows. The total strain energy U stored in an ele- 
ment for any set of arbitrarily varied nodal displace- 
ments u^ and may be written as (refer to Fig. 5 
in the text) 
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xy 
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and 


Consider the continuum element under a state of 
uniform strain e xx (Fig. Al). If Ax represents the 
extension of the element in the x-direction, then. 

Ax/a * e xx ; or, Ax = e xx a. The stresses in the element 
are given by o xx * E^e xx an ^ the force in the 
x-direction in the element is 


+ i v v 4- v 2 *i wv 2 


V v 4 - V 1 + v 3 ~ V 2 )(u 2 


U 1 + u 3 ‘ u 4 ) 



F x = 


a b 
XX 


l e xx b 
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Imposing an equivalent extension in the 
x-direction on the spring element, the force in the 
x-direction at the nodes 2 and 3 are given by K xx Ax 
* K xx e xx a. The total force in the x-direction in the 
element is 



(A10) 

Using Cast i gl i ano 1 s theorem, the nodal forces can 
be obtained as follows: 


F x - 


xl 


+ F 
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Hence from Eqs. (A2) and (A3) K xx = (Enb)/a. 

Similarly, by considering a uniform strain cyy 
in the y-direction, it may be shown that K vy 
- (E 2 2*)/b. 

From Eq. (Al), the stress in the y-direction, 
°y y* 1 n the continuum element due to strain c xx is 
a yy * ^2] E xx s E]2 e xx anc} the force in the 
y-airection is 
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Such derivatives for each joint give the element 
stiffness matrix in the form 


{F} - CKl{d} 


(Al 2) 
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E 1 2 e xx a 
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If K X y is the cross-extensional stiffness coefficient 
of the spring element, then, the force in the 
y-direction due to displacement Ax is: 


where CK1 represents the element stiffness matrix, 
{d} the displacement vector, and {F} the element 
force vector. 

Assembly of the LHR stiffness matrix is done by 
standard direct stiffness methods employed in finite 
element practice. 


F y = 


V x - K xy E xx a 


(A7) 


From Eqs. (A4) and (A5) we have Kxy = El 2 ■ 

Finally, consider the continuum element in a state 

of pure shear (shear strain, e X y), as shown in Fig. A2. 

The shear stress in the element is x = x = E~e . 

xy y x j j xy 

The force in the x-direction on face 1-2 of the ele- 
ment is: 


APPENDIX B - DERIVATION OF THE STRAIN ENERGY RELEASE 
RATE FOR INCIPIENT RUPTURE 

Consider a small virtual increase Ac in crack 
length in an element under a given external load. The 
total potential energy, t r, is given by 

ir = | {u} T nO{u} - {u} T {b} , (Bl) 


F x " T xy a = E 33 e xy a (A8) 

Imposing the corresponding nodal rotations in the 
spring element (Fig. A3), the force in the x-direction 
of the side 1-2 of the element is: 


where {u} is the nodal displacement vector, [K] is 
the element stiffness matrix, and {b} is the nodal 
force vector. 

The energy release rate G is obtained from the 
variation of ir with respect the load, i.e.. 



Hence C = E^ab/2. 

Having obtained the spring stiffnesses in terms of 
the elastic constants, the stiffness matrix of a LHR 


Sir - \ {u} T [Sk]{u} + {Su} T [K]{u} 

- {Su} T {b} - {u} T {Sb} (B2) 

Using equilibrium and the fact that the vector {Sb> 
is null for a given load, the above simplifies to 


6 





I 


Sir - \ (u} T [Sk]{u> (B3) 

Hence, 

G = ' HF= - 1 < U > T Tc CK] < U > (B4) 

if E 6 K] be the change in the stiffness matrix of the 
element due to the incipient fracture, i.e., [ 6 K] 

* [^before - [K] a fter» then G may written in the 
form 



TABLE I. - CONSTITUTIVE MATERIAL PROPERTIES 


Constituent 

Young's 

modulus, 

GPa 

Poisson 1 s 
ratio 

K IC*_ 
MPa m 

SIC fibers 

390 

0.3 

5.0 

Si 3 N 4 matrix 

206 

.3 

4.6 

Interface 

206 

.3 

— 
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FIGURE 3. - LOCAL HETEROGENEOUS REGION (LHR) AND ITS 
COMPONENTS. 
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NO SLIP 



FIGURE 2. - GRAPH DUE TO BUDIANSKY ET AL., FOR CAL- 
CULATING MATRIX CRACKING STRESS. 
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FIGURE 7. - SCHEMATIC OF THE LHR ZONE EMBEDDED IN THE 
FINITE ELEMENT MESH. (NOT TO SCALE). 
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FIGURE 6. - LOCAL FRACTURE MODES IN A LHR ELEMENT. 
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LHR LENGTH, pm 

FIGURE 8. - EFFECT OF LHR LENGTH ON O mal AND <J u|t 

(v f = 0.4, T n = 10 MPa, K M = 4.6 MPaVm). 

1 u IC 
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600 


0 U H , MODEL 



1.00 1.80 2.60 3.40 4.20 5.00 


K M , MPa-s/h 
IC 

FIGURE 10. - VARIATION OF a maf AND O u)t WITH K M 
(T 0 = 10 MPa, v f = 0.4). 
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CRACK OPENING DISPLACEMENT, \xm 



FIGURE 13. - VARIATION OF CRACK OPENING DISPLACEMENT 

u n WITH CRACK LENGTH c (T 0 = 10 MPa, K M = 

0 0 IC 

A. 6 MPaVm, Vj = O.A). 
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EVENT # 

LOAD, MPa 

FAILURE EVENT 

1 

140 

FIRST INTERFACE FAILURE 

2 

150 

INTERFACE FAILURE 

3 

200 

INTERFACE FAILURE 

4 

250 

INTERFACE FAILURE 

4* 

250-295 

VARIOUS INTERFACE FAILURES 

5 

295 

FIRST MATRIX CRACK 

6 

295.5 

MATRIX AND INTERFACE FAILURE 

7 

296 

MATRIX AND INTERFACE 

8 

297 

MATRIX AND INTERFACE 

9 

298 

MATRIX AND INTERFACE 

10 

299 

MATRIX AND INTERFACE 

11 

300-305 

THROUGH MATRIX CRACK 
(BEYOND REGION SHOWN) 

12 

305-500 

EXTENSIVE SECONDARY MATRIX 
CRACKING (NOT SHOWN) 

13 

505 

FIRST FIBER FAILURE 

14 

510 

FIBER 

15 

515 

FIBER 

16 

540 

EXTENSIVE FIBER CRACKING 

17 

550 

ULTIMATE (CATASTROPHIC) 
FAILURE 


(EVENTS MARKED WITH AN ASTERICK (*) DENOTE 
VERTICAL CRACKING WITHIN AN LHR ELEMENT) 

FIGURE 14. - SCHEMATIC REPRESENTATION OF THE PRINCIPAL FAILURE EVENTS IN THE LHR 

ZONE. (T n = 10 MPA, K M =4.6 MPaVm, v f =0.4). 

0 IC f 
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FIGURE 15. - AXIAL LOAD DEFLECTION CURVE FOR VARIOUS 

T n VALUES. (K M = A. 6 MPa\Tm, v f = 0.4). 

0 IC ' 
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1 
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FIGURE A1. - FINITE CONTINUUM ELEMENT UNDER UNIFORM 
STRAIN AND EQUIVALENT LHR ELEMENT. 
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FIGURE A2. - CONTINUUM ELEMENT IN PURE SHEAR. 
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